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Abstract 

Nonoscillatory second order differential equations always admit "special" , principal solutions. 
For a certain type of oscillatory equation principal pairs of solutions were introduced by A. 
Elbert, F. Neuman and J. Vosmansky, Diff. Int. Equations 5 (1992), 945-960. In this paper, 
the notion of principal pair is extended to a wider class of oscillatory equations. Also an 
interesting property of some of the principal pairs is presented that makes the notion of these 
"special" pairs more understandable. 



1 Introduction 



For nonoscillatory equations, the notion of principal solution was introduced by W. Leighton, M. 
Morse and P. Hartman Q. Principal pairs were defined for certain types of oscillatory equations 

in a. 

In this paper, we extend the notion of principal pairs to a wider class of oscillatory equations. 
Under certain conditions, general enough to cover many important special equations (Bessel, Airy, 
etc.) we also study the question of identifying the second member of such a pair when the first 
member is given. We find that among all the solutions, the second member is the one whose zeros 
approximate most closely the zeros of the derivative of the first member of the pair. In these cases, 
this property may also serve as an alternative definition, connecting in a certain sense the Kummer 
and the Priifer transformations and bringing more light to this particular choice of "good" pairs of 
solutions of linear second order oscillatory differential equations. 



2 Basic facts and definitions 



Consider the equation 

y" + q(x)y = (1) 
1 AMS Subject Classification: 34C10, 34C20, 33C10 

Research partially supported by grants from the Natural Sciences and Engineering Research Council (Canada) and 
by Research Grant # 199404 of the Academy of Sciences of the Czech Republic 



1 



on [xo,oo), where q E C°[xo,oo). 

If equation {JxJ) is nonoscillatory, then a solution y\ of (||) is called principal, in accordance 
with the terminology of Leighton, Morse and Hartman (see, e.g., ]5|), if 

y^ 2 (x)dx = oo, (2) 

or, equivalently, if lirn^oo 2/1/2/2 = for every solution y 2 of ([[]) linearly independent of y\. In the 
nonoscillatory case such a principal solutions always exists and is unique up to multiplication by a 
constant factor. 

Now, consider an oscillatory equation (|l]) with two linearly independent solutions 1/1,2/2 nor- 
malized by having their Wronskian 

w = w(y 1 ,y 2 ) = yi{x)y' 2 (x) - y[(x)y 2 {x) (3) 

equal to ±1, i.e., 

\w\ = 1. (4) 

It was proved in Q that if 

Um[y?(x) + y|(x)] = L ^ 0, (5) 



for a pair of solutions y\ , y 2 of an oscillatory equation (|T|) , satisfying (|j) , then the pair y\ , y 2 is 
unique up to orthogonal transformation: 

i.e., if (2/1,2/2) is one such pair then every such pair is of the form (ayi + by 2 , cy\ + dy 2 ) with 
a 2 + c 2 = 1 = b 2 + d 2 , ab + cd = 0. 

In 0, each of these pairs (yi,y 2 ) was called a principal pair of solutions of (jl]). Some sufficient 
conditions for the existence of principal pairs (2/1,1/2) of solutions of (H) with g E C 1 [xo, 00), were 
also given there, based on the fact that v(x) = y\{x) + y 2 {x) is a solution of the corresponding 
Appell equation (see M, ||]) 

v'" + Aq(x)v' + 2q'(x)v = (7) 

on [xo,oo), admitting three linearly independent solutions y 2 (x),y 2 (x) and y\{x)y 2 {x). Here we 
extend the investigation of pairs of solutions to the cases where the L in (||) is or +00. 

We shall need the following important facts concerning the Kummer transformation of the 
second order equations of the form ([!]), derived by O. Boruvka, and summarized in his monograph 
@; see also || Chapter 2]. Consider an equation (|]) (either oscillatory or nonoscillatory) and a 
pair of linearly independent solutions, y\ and y 2 . Denote by w = w(y\, y 2 ) their Wronskian 

w(yi,y 2 ) = yi(x)y' 2 (x) - y' 1 (x)y 2 (x) = w^0. (8) 

Following O. Boruvka, define the (first) phase a of ([[]) with respect to 2/1,2/2 as a continuous function 
a E C°[0, 00) satisfying the relation 

tana = y 1 (x)/y 2 (x), (9) 
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whenever this quantity is defined, i.e., for y 2 (x) 7^ 0- Then a € C 3 [xo,oo), a'(x) / on [xq,oo) 
and 

yi(x) = £\a' {x)\~ l l 2 sina(x), 
2/2(2;) = e|a'(x)| -1 / 2 cosa(x), 



e = ±y— w ■ sign a', 

y{x) = ki\u' {x)\~ l l 2 sin(a(x) + k 2 ) being a general solution of (0) for arbitrary constants k\ and 
k 2 . Moreover, equation (|l[) is oscillatory on [xo, 00) (as x — > 00) if and only if lirn^oo |a(x)| = 00; 
see |]. 

Each solution of (||) is bounded if and only if |a'(x)| _1 is bounded on [xq, 00); see ||. 

Each solution of (||) tends to cis x — > 00 if and. only if lim^ >OQ \a'(x)\ 1 = 0; see §. 

In particular, we have 



fiW + 2/2(2;) 



a'(x) 



on [xo,oo) and for y\ and 2/2 with 10(2/1,2/2) = — 1> we get 

2/i (x) = ±|a'(x)[ -1 / 2 sina(x), 
2/2(2;) = ila^x)) -1 / 2 cosa(x) 



and 



2/ 2 (x) + 2y 2 2 (x) = (a'(x))- 1 . 



3 Principal pairs of solutions — extension 



In the previous section, the number L in the limit relation (|5[) was supposed to satisfy < L < 00. 
This raises the question of extension to the cases L = 0, +00. In both of these cases, the limit 
relation (||) must be supplemented by an additional condition in order to characterize sums of 
squares up to orthogonal transformation. 



Consider a pair yi,y 2 °f linearly independent solutions of ([!]) with Wronskian w(yi,y 2 ) sat- 

1^(^^2)1 = 1- (10) 



isfying 
We can write 

where 
Hence 



y\ = ayi + by 2 , 
V2 = cyi + dy2, 



det 



a b 
c d 



w(yi,y2, 



\ad — bc\ = 1. 



2/i(x) = e|a'(x)| 1//2 [asina(x) + 6cosa(x)], 
2/2(2;) = e|a / (x)| _1//2 [csina(x) + dcos a(x)], 



3 



and 



v ( x ) = y\{x) + yl{x) 

= v(x) (a 2 + (?) sin 2 a(x) + (ab + cd) sin 2a(x) + (b 2 + d 2 ) cos 2 a(x) 



(11) 



In case v(x) — > L with < L < +00, it is clear that v(x) approaches a finite limit if and only if 
a 2 + c 2 = b 2 + d 2 = K > 0, and ab + cd = 0. However, in case L = the condition u(x) — » 
implies that U(x) — > 0, regardless of the values of a, 6, c, d. In case L = 00, v(x) — > 00 implies that 
u(x) — > 00, for any values of a, 6, c, d satisfying a 2 + c 2 = b 2 + d 2 > \ab + cd|. Thus the condition 
(|5|), with L = or L = 00, does not characterize the pair 2/1,2/2 U P to orthogonal transformation. 
In as much as this excludes some important examples, such as the Airy equation (q(x) = x) and 
the Cauchy-Euler equation (q(x) = ^y 2 /x 2 ), we now present a theorem which characterizes the pair 
2/1)2/2 by a condition on v': 

Theorem 1. Let (|7p be an oscillatory equation and let 2/1 , 2/2 be a pair of linearly independent 
solutions normalized by the unit Wronskian \ w(yi, 2/2)! = 1- Let v{x) = y 2 {x) +2/|( x ) an d suppose 
that 

lim v'(x) = K, < K < 00. (12) 

x — >oo 

Then the pair 2/1,2/2 is unique up to orthogonal transformation. 

Proof. Let a denote a phase of equation (||) corresponding to 2/1 and 2/2- In accordance with 
§2, we have a G C 3 [xq, 00), a'(x) / on [xq, 00), and 



lim \a(x) 

x— >oo 



00, 



(13) 



because ([[]) is an oscillatory equation. Moreover, 

2/1 (x) = e|c/(x)| -1 / 2 sina(x), 

2/2(3?) = e|a'(x)| -1 / 2 cosa(x), 
e = ±1. Hence v(x) = y\{x) + 2/I( x ) = l a/ ( x )l _1 - Now, for a pair yi,y 2 , we get, from (11) 



^(x) 



V X 



(a 2 + c 2 ) sin 2 a(x) + 2(ab + cd) sin a(x) cos a(x) + (6 2 + d 2 ) cos 2 a(x) 



+ ■ 



1 



a'(x)\ 



i.e., 



(a 2 + c 2 - 6 2 - d 2 ) sin 2a(x) + 2(a6 + cd) cos 2a(x) 



U'(x) =v'(x)F(x) + G(x) 



a x 



(14) 



where 



F(x) = (a 2 + c 2 ) sin 2 a(x) + 2(ab + cd) sina(x) cos a(x) + (b 2 + <i 2 ) cos 2 a(x) 
is a bounded function on [xo, 00), and 

G(x) = K\ sin 2a(x) + cos 2a(x) , 

where K\ = a 2 + c 2 - b 2 — d 2 , K 2 = 2{ab + cd). Now we show that if 

= K, \K\ < 00. 

i.e., if (|12|) holds for yi,y 2 , then iTi = = 0. 



lim t; (x) 

x— >oo 



(15) 
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We have 

2F(x) = {a 2 + b 2 + c 2 + d 2 ) - K x cos 2a{x) + K 2 sin 2a{x) . 
Then, using (111), we get 

2v'( x ) = v'{x)[a 2 +b 2 + c 2 + d 2 ] + cos 2a(x)[K 2 - v'(x)Ki] + sin2a(x)[K 1 + v'(x)K 2 ]. (16) 

Letting x — > oo, we get i^i = K 2 = 0. Hence a 2 + c 2 = 6 2 + d 2 = K% ^ 0, and ab + cd = 0, so 

/a iWa c\ /X 3 \ 
[ c d J ' [b d ) [ K 3 J' 

From the normalization of the Wronskian, we get \ad — bc\ = K3 = 1; hence (yi,y 2 ) is unique up 
to orthogonal transformation. This completes the proof of Theorem 1. 

3.1 An extended definition 

Definition. The pair yi,y 2 will be called principal for an oscillatory equation (fjj if v(x) := 
y 2 (x) + y 2 (x) satisfies either 

lim v(x) = L, (0 < L < 00), (17) 

x — >oo 

or 

lim v'(x) = K, (0 < K < 00). (18) 

x — >oo 

Due to Q, and Theorem 1, this pair is unique up to orthogonal transformation. 

Since all principal pairs for an oscillatory equation (if they exist) differ only by an orthogonal 
transformation, the expression 

v ( x ) = yi( x ) + 2/10*0 

remains the same for all principal pairs. 

The following sufficient condition for the existence of principal pairs is simply a restatement 
of a result of P. Hartman || Theorem 20. lo]: 

Corollary 1. Let q G C 1 [xo,oo), q'(x) > 0, q"(x) < and lim^oo q(x) = 00. Then ^) has 
a principal pair of solutions y\ , y 2 such that 

v(x) > 0, lim v(x) = 0, 

v (x) < 0, lim v'(x) = 0, 
lim v"(x) = 

x— >oo 

and 

lim v'"(x) = 0. 

x^oo 



Another sufficient condition for the existence of principal pairs follows from another result of 
Hartman || Theorem 22. Iq]: 
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Corollary 2. Let q G C 2 [x ,oo), q'{x) < 0, q{x)q"(x) — 3q'(x) 2 > 0. Then ^) has a principal 
pair of solutions 2/1,2/2 such that 

v(x) > 0, 



lim 

x— >oo 



d 



<l \ x ) [fa] v ( x ) 



0, n = 1,2,3, 

and v approaches a finite limit or oo according as q approaches a positive limit or 0. 



3.2 Examples 

Example 1. The generalized Airy equation 

This is an example with L = 0. The equation in question is 

y" + (2u)- 2 x l ^- 2 y = 0, (19) 

where we suppose that < v < 1/2. The usual Airy equation corresponds to v = 1/3. With 
the usual notation for the Bessel functions M, the pair x 1 / 2 J v (2vx 1 /^), x l / 2 Y v (2vx 1 ^ 2u " > ) satisfy 



(|l|). In view of the known result |, p. 446] that, for < fi < 1/2, the function t[J 2 (t) + Y 2 (t)] 



increases to 2/tt on (0,oo), we see that 

v(x) = x[J 2 (2ux 1 /^) + Y 2 {2ux 1/{2v) )\ 

approaches as x — > oo. In fact, for \ < v < \, v(x) is completely monotonic (i.e., its successive 
derivatives alternate in sign) Q Theorem 5.1 (i)]. The asymptotic formula 



x[J 2 (2ux 1 ^) + Y 2 (2vx 1 '^)] = —x 1 - 1 '^ [l + 0(x~ l H , (20) 



& p. 224 (5)] shows that v(x) and v'(x) approach as x — > oo. Thus, by Theorem 1, the pair 
x l / 2 J v {2vx l /^),x 1 / 2 Y u (2ux 1 /^) constitute a principal pair for (||). 

Example 2. 

The equation 

y" + x- 1 y = 0, (21) 

has a principal pair x 1//2 J\(2x 1 ^ 2 ), x 1//2 Y"i(2x 1 / 2 ). Here v(x) — > oo, but v'(x) — > as x — > oo, so the 
hypotheses of Theorem 1 are satisfied again. 

Example 3. The Cauchy-Euler equation 

This is an example with L = oo, to which Theorem 1 is applicable. The Cauchy-Euler 
equation is 

y" + \v = o. (22) 

x A 

In case j 2 > 1/4, the equation is oscillatory and a principal pair of solutions is given by yi(x) = 
x 1 / 2 sin(slogx), y2(x) = x 1 / 2 cos(s log x), where s = ^j 2 — 1/4, so v(x) = x and v'(x) = 1, so the 
hypotheses of Theorem 1 are satisfied. 

It is clear that every pair 

Vi = ay i + by 2 , 
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y 2 = cyi + dy 2 , 

with a 2 + c 2 = b 2 + d 2 > \ab + cd\ has the property that v{x) — > oo as x — > oo, but it is only in the 
case a 2 + c 2 = c 2 + (i 2 , a& + cd = that u'(x) approaches a finite limit. 



4 Properties of principal pairs 



Suppose that an oscillatory equation (1) admits a pair y\,y 2 , for which 

lim v'(x) = lim (yf(x) + 2/ 2 (x))' = 0. 



(23) 



Due to the extended definition in §3.1, the pair 2/1,2/2 is a principal pair. Let a denote the phase 
of (1) corresponding to 2/1,2/2- Then 

2/1 (x) = e|a'(z)| -1 ' 2 sin a (a;), 

2/2^) = ela^x)!^ 1 / 2 cosa(x), 
e = ±1, 

a £ C 3 [xo, 00), c/(x) 7^ on [xo, 00), 

= yi( x ) + yiOO = la'C^OP 1 - 



Condition (G3h gives 



Consider the first derivative of 2/1 : 



lim 

x^oo (a'(x)) 2 



.I| a '( x )|-i/2^M sina(x) + \ a '(xp 2 



2 a'lxj 
All of its zeros are at points Xj where 

cos a(xj) 1 a"{x,j) 



cos a [x ■ sign a 



sina(xj) 2(a'(xj)) 2 



From (l24|) we have 



(24) 



(25) 



\a(xj) + 7r/2 — j7r sign a'| — » 0, as j — > 00. 
A general solution of (||) can be written in the form 

y(x; ci, c 2 ) = ci\a' (x)\~ l l 2 sin(a(x) + c 2 ). 

Its zeros are at points Xj, where 

a{xj) + c 2 = jir sign a', j = j , j + 1, . . . . 

Hence 

y(x; ci, it/2) = ci\a' (x)\~ l l 2 cos a(x) 

is that solution (up to constant multiple) whose zeros are approached by the zeros of the derivative 
of 2/1. However, after the normalization c\ = ±1, y\ together with y(x, ±1, 7r/2) = ±2/ 2 (x) represent 
a principal pair of (0). 
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We may summarize our considerations in the following Theorem: 

Theorem 2. If an oscillatory equation ^) has a pair of solutions 1/1,2/2 for which 

lim v'(x) = 0, (26) 

x— >oo 

then they have the property that the difference of the zeros of the derivative of one of them and the 
zeros of the other tend to as x — > 00. In other words, among all sequences of zeros of solutions 
of the best approximation to the zeros of y[ is given by zeros of y 2 - 

Remark. Evidently, the conditions in Corollary 1 are sufficient to guarantee the conclusion 
of Theorem 2. Further sufficient conditions can be obtained by using Hartman |], Theorem 18. lo]. 
In fact, this Theorem, together with his Theorem 20. lo, show that the requirements 

q' > 0, q" < 0, and lim q(x) = q(oo), where < q(oo) < 00 

x — >oo 

are sufficient to imply the conclusion of our Theorem 2. 

The property described in Theorem 2 links two transformations of second order equations. 

The first is the Kummer (sometimes called Stackel) transformation consisting in the change 
of independent and dependent variables of the solution z(t) into the solution y{x) expressible in 
the form 

y(x) = f(x)z(h(x)) 

or, better, 

Vi( x ) = f{ x )z\{h{x)) 

and 

2/20*0 = f(x)z 2 (h(x)), 

for two linearly independent solutions y±, 1/2 and z±, Z2 of the transformed equations. In particular, 

yi(x) = f{x) sin h{x) 

and 

y 2 (x) = f(x)cosh(x), 

when the equation z" + z = is taken as canonical (as can always be done). The second is the 
Priifer transformation, considering again a two dimensional underlying space for two functions, now 

yi(x) = y(x) = p{x) s\mp(x) 

and 

2/2(2?) = y'(x) = p(x) cos ip(x). 
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